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Abstract. The indefinite integral of the homogenized Ornstein-Uhlenbeck process is a 
well-known model for physical Brownian motion, modelling the behaviour of an object 
subject to random impulses [L. S. Ornstein, G. E. Uhlenbeck: On the theory of Brownian 
Motion. In: Physical Review. 36,1930, 823-841]. One can scale these models by changing 
the mass of the particle and in the small mass limit one has almost sure uniform conver- 
gence in distribution to the standard idealized model of mathematical Brownian motion. 
This provides one well known way of realising the Wiener process. However, this result 
is less robust than it would appear and important generic functionals of the trajectories 
of the physical Brownian motion do not necessarily converge to the same functionals of 
Brownian motion when one takes the small mass limit. In presence of a magnetic field the 
area process associated to the physical process converges - but not to Levy's stochastic 
area. As this area is felt generically in settings where the particle interacts through force 
fields in a nonlinear way, the remark is physically significant and indicates that classical 
Brownian motion, with its usual stochastic calculus, is not an appropriate model for the 
limiting behaviour. 

We compute explicitly the area correction term and establish convergence, in the small 
mass limit, of the physical Brownian motion in the rough path sense. The small mass limit 
for the motion of a charged particle in the presence of a magnetic field is, in distribution, 
an easily calculable, but "non-canonical" rough path lift of Brownian motion. Viewing the 
trajectory of a charged Brownian particle with small mass as a rough path is informative 
and allows one to retain information that would be lost if one only considered it as a 
classical trajectory. We comment on the importance of this point of view. 



1. Introduction 

Newton's second law for a particle in with mass m, and position x = x (t) , (for 
simplicity: constant) frictions 01,02,03 > in the coordinate axis, subject to a (3- 
dimensional) white noise in time, ^, where ^ = .^(t) is the distributional derivative of 
W, a (mathematical) Brownian motion or Wiener process, reads 

(1.1) mx = —Ax + ^ 
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where A = diag (ai, 02, 03)- Orthonormal change of coordinates imphes that the "correct" 
assumption for A is to be symmetric with strictly positive spectrum, 

a (A) C (0,00) . 

The process x (t) describes what is known as physical Brownian motion. It is weh-known 
that in smah mass regime, m « 1, of obvious physcial relevance when dealing with 
particles, a good approximation is given by (mathematical) Brownian motion; to see this 
formally, it suffices to take m = in (II. ID in which case Ax is a standard Brownian motion 
in M^. 

Let us now assume that our particle (with position x and momentum mx) carries an 
electric charge q ^ and moves in a magnetic field B which we assume to be constant. Re- 
call that such a particle experiences a sideways force (" Lorentz force" ) that is proportional 
to the strength of the magnetic field, the component of the velocity that is perpendicular 
to the magnetic field and the charge of the particle. 

FLorentz = qx X M. 

When B is constant, which we assume for simplicity, the Lorentz force experienced by the 
particle (at time t) can be written as as linear function of velocity x = x (t), namely qBx 
for some anti-symmetric matrix B. In other words, the dynamics for physical Brownian 
motion in a magnetic field are given by 

mx = —Ax + qBx + ^ 
= -Mx + C) 

where M = (A + qB) is such that all eigenvalues of M G M"^" have strictly positive real 
part (one may still think n = 3, but the subsequent analysis works for any dimension 
n). Note that these second order dynamics can be rewritten as evolution equation for the 
momentum 

p (t) = mx (t) , 

indeed, 

p = -Mx + C = -—Mp + W, 
m 

and we shall take this point of view when rewriting the dynamics in term of standard 
stochastic differential equations. As before, when m = 0, the process Mx = is a bona 
fide (i.e. mathematical) 3-dimensional Brownian motion and one may think that little has 
changed, appart from the covariance matrix of the resulting Brownian motion. And indeed, 
writing £, = W, and assuming xq = Wq = 0, it is easy to see that 

Wt — Mxt = / {—Mxs + £,s)ds = / p = {pt — Po) — > as mass m — 0. 
Jo Jo 

Note that 

Po = mx (0) — )• 

as m — 7- 0, whenever initial velocity remains uniformly bounded, so the statement here is 
that — 7- as m — 7- and one can easily see that this convergence is uniform (we are 
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only interested in a fixed time horizon, say [0, T]). However, the momentum may have 
quite non-trivial effects as control. Recalling that controlled system (differential equations, 
integrals ...) behave in a robust fashion precisly under rough path metrics, the essence 
of which has been briefly summarized for the reader's convenience in the appendix, the 
following lemma tells us that momentum does not at all converge to a trivial limit. The 
situation is reminiscent a well-known deterministic example, taken from [18], where the 
path 

Zt := exp (27rztm) G C ~ 



m 



converges to a non-trivial "pure area" rough path as m — t- 0. (See, however, remark [T3l 
where we emphasize the special role of Brownian motion, notably the need for its intrinisic 
irregularity, in the proposition and theorem below.) 

Proposition 1. One has (the iterated integrals are understood in ltd sense) 

{Ps,u ^s,t) ■■=[ dp, / dp(E)dp] ^ {t-s)] as m ^ 0. 



and C is the (symmetric) n x n matrix defined by 

C= / e-'^^'e-^'^'ds. 





More precisely, the convergence holds in the following strong sense: for any a G (1/3, 1/2), 

\PsA \¥s,t-MC^it-s)\ 
1+ \l \ 

s,te[0,T] F - si s,iG[0,T] F - s| 

(This is precisely what is meant by convergence in a-Holder rough path metric.) 



(1.2) sup j2 — ^^-|^ -|- sup - — i — ^ — >• m L'^, as m ^0. 



In view of this proposition one suspects (correctly) that the convergence of physical 
Brownian motion to Brownian motion is also non-trivial if one thinks how Brownian par- 
ticles act as controls, i.e. as rough paths. More specifically, one expects a limit in which 
Levy's stochastic area is perturbed by a term proportional to 

{MC -CM*)/2, 

the (anti-symmetric) matrix which effectively described the pure area (rough path) limit 
of the previous proposition. Let us insist, however, that such a statement (i.e. the content 
of the following theorem) is not a corollary of the above since, in general, 

y (Z - Z") d (Z - Z") ^ ^ j Z'^dZ'^ j ZdZ. 

Indeed, if one thinks of Z", Z as rough paths, their (formal) iterated integral are meaningful 
by definition of a rough path. In contrast, the iterated integrals of Z agains Z" will not 
even be, in general, well-defined. We are now ready to state our main result. 
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Theorem 1. Let M be a square matrix in dimension n such that all its eigenvalues 
have strictly positive real part. Let W he a n-dimensional standard Brownian motion, 
m ("mass") as strictly positive scalar and consider the stochastic differential equations 

(1.3) dX = — Pdt (position) 

m 

(1.4) dP = — ^-MPdt + dW (momentum). 

m 

with (for simplicity only) zero initial position and momentum. Let W = (W^W) he the 
natural rough path lift ofW , where Ws,t = Jl (Wr — Ws)<SSiodWr. Define also W = (W,W), 
where 

(1.5) W,,t = Ws,t + {t- s)^{MC -CM*), 

and C is as in the previous proposition. Then, as m ^ 0, MX converges to W in L'' and 
a-Holder rough path topology, for any q > 1 and a £ (1/3, 1/2). More precisely, we have 



(l.bj sup i i h sup 



Jl MXs,r ® d{MX)r - W,,t 



in L'^ 



as m — )• and this convergence is of rate arhitrarily close to 1/2 — a. 



Remark 1. Ln view of the tensorial transformation hehaviour of iterated integrals, ^.6\) 
is plainly equivalent to 



Xs,t ^ M-^Ws,u 
* 1 



Xs,u®dXu ^ M-^Ms,t{M~^)* + {t-s)-{C{M-^Y -M-^C); 



in a-Holder rough path metric. 
Remark 2. One has 

W = W 

if and only if M is symmetric. To see this, note that MC is symmetric (hence equal to 
CM*) if and only if M is symmetric, using symmetry and invertihility of C . 

Remark 3. The framework of Gaussian rough paths [7] and [9] Ch. 15], which plays keyrole 
in non-Markovian Hormander theory [3l |3] and some recent hreakthroughs in non-linear 
SPDE theory ITU] is only applicable if M is symmetric since then one can diagonalize 
the dynamics such as to have Gaussian driving signals with independent components. Ln 
this case, upon checking some uniform variation hounds on the covariance, it could he used 
to see the convergence il.6\) to W, the "natural" rough path lift of Brownian motion; [Hi 
Ch. 15] . But since W is not natural when Anti (M) ^ 0, we here also provide an explicit 
example which illustrates the necessity of the assumptions put forward m [71 19] . 



PHYSICAL BROWNIAN MOTION IN MAGNETIC FIELD AS ROUGH PATH 



5 



Remark 4. If M is normal, i.e. 

[M, M*] = 0, 

then the difference between W and W has a somewhat simpler expression. Indeed , we 
compute 

C = ^Sym(M)-\ 
and since C commutes with M , we get 

(1.7) \{MC - CM*) = i Anti(M) Sym(M)-^ 
Similarly, the area correction for X is seen to simpify to 

(1.8) ]^{C{M-^)* - M~^C) = Anti(M-i) Sym(M)-\ 

Remark 5. It would he possible to "unit" the above proposition and theorem in saying that 
the physcial Brownian motion is a "good" approximation in the sense of [6j to the limiting 
rough path W. 

Remark 6. Incidentally, the rate "anything less than 1/2 — a" is also the rate of conver- 
gence (in a-Holder rough path metric) for piecewise linear "Wong-Zakai" approximations 
to Brownian motion obtained in [ll] and optimally of this rate is known. In both cases, 
these rates are obtained as a rather mechanical consequence of good moment estimates, cf. 
[Hi Thm A. 13], so that we also believe our rates to be optimal. 

The following example is taken from Pavliotis-Stuart, Hairer \22\ Section 11.7.7]. We 
note that the PDE argument (based on multiscale techniques) offered in [22[ Section 11.7.7] 
can only give convergence of the finite-dimensional distributions, tightness - especially in 
rough path topology - will require additional and non-trivial estimates (which are implied 
for our work below). And of course, strong convergence, available to us because of a natural 
coupling of physical Brownian motion X and W = is out of reach with PDE methods. 

Example 1. Take n = 2, a G R and 

M = Id-a(^^^ 

so that 

l + a^ \ a 1 J 

Note Sym(M) = Id. Then, from il.8\) . the area correction of X , in the m — t- limit, 
equals 

where 61,62 denotes the standard basis o/R^. This agrees precisely with [22, (11.7.28c)]. 
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Let US conclude this introduction by noting that the example of physical Brownian 
motion under influence of a magnetic field as a rough path has some history. Indeed, it 
appeared as motivation (but without much details) in several presentations, including some 
by the last-named author in early 2000. 

2. Proof of theorem [T] 

We first give a proof of the rough path convergence without rates, based on the ergodic 
theorem. The adaptions which lead to the announced rates (and also bypass the use of the 
ergodic theorem!) are then explained in details in the remark following the proof. 

Proof. In order to exploit Brownian scaling, it is convenient to set 

m = 

and then as rescaled momentum, 

= P/e. 

We shall also write = X, to emphasize dependence on e. We then have 

(2.1) dY^ = -e-'^MY^dt + e~^dW 

(2.2) dX"" = e'^Y'dt, 

By assumption, there exists A > s.t. the real part of every eigenvalue of M is (strictly) 
bigger than A. For later reference, we note that this implies the estimate 

I exp(-rM)| = 0(exp(-AT)) 

as r — )• oo. For fixed e, define the Brownian motion 

W. = eW,-2., 

and consider the SDKs 

dY = -MYdt + dW 
dX = Ydt; 

note that the law of the solutions does not depend on e. Furthermore, when solved with 
identical initial data as (12.ip - (l2.2j) . we have pathwise equality 

(2.3) {Y^X'/e) = Z-2.,X,-2. 

Thanks to our assumption on M, Y is ergodic; the stationary solution has (zero mean, 
Gaussian) law 

i/~AA(0,C) 

for some covariance matrix C. To compute it, write down the stationary solution 

y^stat ^ e-^'^'^'-'^^dWs] 
J — oo 



Since 
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for each t (and in particular t = 0) the law of Y^^^^ is precisely i^. We then see that 

/"O poo 
oo Jo 

sup E|y/| < oo, 

0<t<oo 

it is clear that 

eY,-2t = eVf ^ 
in uniformly in t (and hence in L'^ for any q < oo). Noting that 

MX^ = Wt-eYo% 

the first part of the proposition is now obvious. Moreover, by the ergodic theoren{3, 

(2.4) 1^ f{Yf)dt ^ 1 1 fivHdy), in L\ any q < oo, 

for all reasonable test function /. We shall only use it for quadratic^. Using 

dX^ = e-^Y^dt 

we can then write 

ft 



' s 

(■t 



[ MXl d{MX^)s = [ MXl ®dWs-e [ MXl ® dY^ 
Jo Jo Jo 

= [ MXI(^dWs-MXf(^{eYf) + e I d{MX^)s(^Y^ 
Jo Jo 

= [ MXl0dWs-MXf0{eYf)+ [ MY^ (^Y^'ds 
Jo Jo 

Ws(^dWs-Q + t I My0 yv{dy) 







t 

Ws dWs + tMC 



yo,t + t ( MC - ll 



where the convergence is in L'^ for any q > 2. By considering the symmetric part of the 
above equation, 

^(MXf) ® (MXf) ^ ^Wt (^Wt + Sym (^MC - , 



^See e.g. [23] (p.421) or [13 (p.409) 



^The ergodic theorem in the references we have cited only apphes to bounded /, but it is easy to see 



by a truncation argument that (|2.4p still holds as long as (/(Vs))s>o is bounded in any L'' and u has finite 
moments of all order. 
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we see that 



MC-\l 



has symmetric part 0, i.e. is antisymmetric, and hence also equals 

^ (MC - CM*) . 
This settles pointwise convergence, in the sense that 



MXf, I MXl®d{MX^ 





In view of |9j. Proposition A. 15, the rough path convergence ()1.6p will follow once we have 
checked that for any g < oo, 



sup E[\Xlt\'^] < Cq\t-S 
eG{0,l] 



sup E 

£G(0,1] 



Xt . ® dX' 



< Cq\t-S\'^. 



First, since X is Gaussian, it follows from integrability properties of the first two Wiener-Ito 
chaos, that it is enough to show it for g = 2. Secondly, we note that the desired estimates 
are a consequence of the following : 



(2.5) 
(2.6) 



E 



X 



s,t 



E 



X^.(R)dX 



< (const. )|f — s|, 

< (const.) \t — , 



where the constants must be uniform over t, s G (0,oo). Indeed, this follows directly from 
writing 



E 



\^s,t\ 



E 



eX,- 



< (const. )e \e t — e s = (const. )|t — s| 



(note uniformity in e), and similarly for the second moment of the iterated integral. 

In order to check (|2.5p . it is enough to note MXs^t = Wg^t — Ys,t, combined with the 
estimate 



E 



s,t\ 



E 



-M(t-s) 



+ / rr(e 



-Mu-M*u 



)du 



< C{M)\t-s\, 
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E 



< 



E 



t fU 



YlYUrdu 



-^r u q v 



'^{r<u;q<v}dr du dq dv 



E 



r ^ u 



E 



Y^Y^ 

q ^ V 



+ 



E 



Y^Y^ 

1 J. 1 g 



E 



Y^Y^ 



r V 



E 



u q 



dr du dq dv 



Yr ® Yu 
Yr ® Yu 



dr du 



l|r<^}dr du 



where we have used the fact that Y is Gaussian ("Wick's formula"). But for r < u, 

E \Yu\Yr] = e-^("-'')i;, so that 



E 



E 



^{r<u}dr du 



Yr(S)Yu 

Yr e-*^^"-'^)!;. 



^{r<u}dr du 



s \Jr 



Yr 



dr] <C3(^-s)^ 



recahing that |exp(— rM)| = 0(exp(— At)), and (j2.6p is proved. 



□ 



Remark 7. (Rates) The use of the ergodic theorem can be avoided by checking "directly' 
that 



Yf Yfds ^tC = tEi Fq'*"* ^ I'o'*"* 



in L'i Vg < 00, with a proof similar to the proof of the inequality ()2.6p . Assume for 
simplicity that Y^ = Ygj^/i is started with Yq ~ v; note then that E (g) y/ = tC . 
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Furthermore : 

r-t 



E 







[0,t] = 



[0,t]2 



E [y/'V/'-^y/'^y/'-'] ds dr 

E [y/- V/'^] E [y/'^y/'^] ds dr 



+2 [ E [y/'*y/'^'l E [y/'V/'^'l ds dr 

< -t^ + f (e [y/'V/'^] ^ + E [y/'^y/-^] ^) ds dr. 

< Cf£ + 4 [ E [y/ ® y/]2 1 1,<^| ds dr. 



A'^ow noie i/iai for s < r, 
so that 

and we finally obtain 
ft 



E[y/|y/] =e-' 'M(r-s)ye^ 

E[y/®y/]2 = o(e-!^ 



E 







Y^'''Y,^'^ds - tCij 



E 



t N 2- 

Y!''Y!'^ds 







t / ;>oo 



'^^dr] ds = (const. Je"^. 



< (const.) 

We have thus proved that the L'^ -convergence (and by Gaussian properties L'^ -convergence 
for any q < oo) JqY^ y/ds — )• tC is in fact of order e. Actually, from here on it is not 
difficult to establish convergence rates of il.6\) . From [9l Thm A. 13] and the work already 
done in the previous proof (reduction from g*^ moments to second moments is immediate 
by Gaussian chaos) it will suffice to check 

|2" 



(2.7) 
(2.8) 



E 



\Ws,t-MX'^ 



E 



s,t - / MXl, d{MX')r 



< {const.)5{e)\t- s]'^, 

< {const.)6 (e) \t - sl^/^ , 



for fixed 1/3 < (3 < 1/2 and 6 = 6{£;f3). As long as S {e) ^ with e, this is also the 
convergence rate in M.(?) . for any 1/3 < a < /?. A short computation shows that 



E 



Ws t - MX 



s,t\ 



E 



MXl^. d{MX^ 



< C\t-s\^/^ and <Ce, 



< C\t — s\ and < Ce. 
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By (geometric) interpolation, with exponent 2(3 < 1 and 1 — 2/3 

1 

|2 



E 



W, t - MX 



s,t\ 



E 



>s,t 



< Ce^-^P\t-s\'r 



and we obtain (and then L'^ , any q < oo) convergence rate 6 (e) = e^^"^^ . In other 
words, given a G (1/3,1/2) we have rate arbitrarily close to 1 — 2a; by a Borel-Cantelli 
argument this is also the almost-sure rate, say, along e = 1/n. Note that by working in 
the stronger topology (a f 1/2/, one loses on the rate of convergence. Also, since "level-2" 
rough path theory imposes a > 1/3, an upper bound for the best possible rate (in a rough 
path metric!) is given by 1 — 2/3 = 1/3. (The situation is very similar to the rate of 
convergence of piecewise linear approximations to Brownian motion (as rough path), see 

m. 



3. Applications and remarks 



We conclude this note with a number of applications and remarks. 

Remark 8. Write the (anti- symmetric) matrix MC — CM* as Yli<j1i,j[^i'^j]- ^ 
consequence of our main theorem, we have the following convergence result for ODEs driven 
by "physical" Brownian motion, in the zero mass limit. Let Vq G Lip^, V = {Vi, . . . , Vn) G 
]^-p2+(5 j^j, some 5 > 0, be vector fields on M^. Let be the solution to the SDE (actually, 
random ODE) 

d 

dY^ = VQ{Y^)dt + Vi{Y^^)dX^'^\ Y^ = xe M^ 

i=l 

Then as e — )• 0, Y'^ converges to the solution Y of the following ltd SDE 

n 

dY = Vq {Y) dt + Y, Vi{S°)dBi, 
1=1 

where {Vi, . . . ,Vn) = {Vi, . . . , and 

i=l l<i<i<" 

Indeed, given our main result, this is a simple consequence of rough path stability and the 
impact of higher order perturbations on RDEs, see e.g. ^ Ch. 12], combined with the 
usual Ito-Stratonovich correction. 

Remark 9 (Magnetohydrodynamics). A (physical) system - as described by the above 
differential equations - which is driven/controlled by a single, charged Brownian particle in 
a magnetic field may not appear to be a standard situation in applied science. However, it 
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is not hard to think of a system being influenced by a cloud of such particles. If N such 
particles move independently, our main theorem applies immediately in dimension n = 3N. 
Of course, for N » 1 one needs to incorporate interactions between the particles. In fact, 
the movement of a cloud of charged particles will effect very much the magnetic field itself. 
In our model, the matrix M should then be allowed to depend on the (bulk) behaviour of 
the N particles. Much more work will be necessary to give a proper rough path analysis 
of this situation, we do remark however that first mean field results for rough differential 
equations have been obtained by Cass-Lyons [5] . It is quite conceivable that our rough path 
perspective then provides a very novel point of view for magnetohydrodynamics. 

Remark 10. Similar correction terms appear when one considers limits of one-forms inte- 
grated against . The rough path correction matters, for instance, if one wants to make a 
change of measure and represent the stochastic integral in the Girsanov factor in terms of 
, rather than in terms of the underlying Brownians. (This type of representation plays 
an important role in "robust" filtering, path sampling of conditioned diffusions \i2\ 111] and 
related topics.) 

Remark 11. The signature of a path is the sequence of its iterated integrals against itself. 
For (deterministic) paths of bounded variation, it fully characterizes the path up to "tree- 
like" equivalence [13j. In a similar spirit, the expected signature of a processes (in the 
sense below) characterizes the essence of its law, at least when it comes the solution of 
stochastic differential equations; the so-called cubature method is based on this [20]. In [15] 
expected signatures of many stochastic processes are considered. By either specializing these 
considerations or a direct computation one can see that, as e — )• 0, the expected signature 
of MX^ converges to 



Remark 12. In this paper we have only considered the case where M was constant. It is 
however natural to consider the case where the dynamics (i.e. friction, magnetic field and 
covariance matrix of the brownian term) depend on the position X of the particle. This 
leads to consider the coupled system of SDEs 



It is then not difficult to show that the paths X^ — >• X pointwise, where X is solution 



dX^ = [D^{QM-^) : C] {X^)dt + {QM-^){X^)dW, 
where [D^{QM-^) : C]' := Y.k,idk{QM-^)\C^^ and C = C{x) is still defined by the same 




dX^ 



dY' 



-£~^M{X')Y'dt + e'^dWt, 
£-^Q{X^)Y^dt. 



to 



formula. 
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As for the iterated integral, we obtain a similar correction except that it is now state 
dependent (and thus random): 

f X^®dX^^ I X^C^odX^ + \ [ Q[C{M-^Y - M-^C)Q*{X^,)ds. 
Jo Jo ^ Jo 

The computation of these pointwise convergences is close to the beginning of the proof of 

Theorem 1, but instead of using the ergodic theorem one should notice that for small e 

ft ft 



f{XlY!)ds ^ l^ij f{Xly)uc{Xldy))ds, 



as the process evolves at a much faster time-scale than X^ (here ~ should mean that 
the difference is small in Lp'-norm). The detailed verification of convergence in rough path 
sense is technical (mainly, because one looses the Gaussian and Markovian structure of 
y^) and left to subsequent work. 

Remark 13 (On the role of Brownian roughness). Let us return to but now with 

Brownian motion replaced by a (deterministic) path 7 defined on [0,T]. That is, we consider 
the evolution 

mx = —Mx + 7. 

Of course, even when 7 fails to be differentiable this equation is well-defined in the distri- 
butional sense, thanks to the additive structure of the noise. As before we assume that M 
has an antisymmetric part which therefore triggers rotation and thus effects the area. One 
may wonder how "rough" the driving force (now assumed deterministic) needs to be to see 
some non-trivial area correction in the limit. As we now show, the roughness of Brownian 
motion -with (almost) 1/2-Hdlder sample paths - is crucial and no area corrections arises 
when 7 is Holder with exponent greater than 1 /2 . 

Proposition 2. Assume 7 is a-Holder. Then, as m — t- 0, Mxq^. converges to 79 . in 
j3-H6lder topology for any (3 < a. In particular, when a > 1/2 it follows that 

Mxo^t <X) d (Mxt) / . (g) dj as m —?- 
' Jo ' 

where the integral on the left-hand- side is understood in Young sense, convergence may be 

understood uniformly on [0,T] (and actually in (3-Hdlder rough path sense). 

Proof. By interpolation, it is enough to establish pointwise convergence of Mxq,- to 79 . in 
conjunction with uniform Q-Holder bounds. Equivalently, we want to show that, pointwise 
and with uniform a-Holder bounds, 

z := 7o^. — Mxq^. — )• as m — )• 0. 

Note that 

M , 

dzf = ztdt + a7f 

m 

from which we see, writing z^- t = zt — Zg as usual, 

t 



s 
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The last integral is a Young (actually Riemann-Stieltjes) integral, for its integrand has 
finite variation. To see this, note that for s <t, |e~*^ — e~*^| < Ce~^^ {{t — s) A 1) where 
A is, as in previous sections, a lower bound on the real part of the spectrum of M. We may 
then estimate, for any subdivision 0<so<...<SAr<oo 



E 



i 

< e"^" E Hn<s,<n+1} ((^i+l - Si) A 1) 
n=0 i 

oo 

< 2(7^6"^" < 



In particular, it follows that 



sup 

0<m<l 



n=0 



_ M 
e m ' 



l-var;[0,T] 



OO. 



< OO. 



We now address pointwise convergence. Since 2;o = we can estimate, whenever Q < 5 < 
t < T, 



f 

Jo 



< e 



pt-S 

Jo 



_M_ 

e m 



iit-S)-r}dj 



+ 



f 

Jt-t 



< Ce 



7lla-H61;[0,T] 



T" 1 + 



+C||7lU-H6I;[0,T]'5" 1 + 



< C (e-^^ + 5") , 



e ™ 



l-var;[0,T] 
l-var;[0,T] 



where C is a constant which does not depend on m. Taking 

am 



Alog(l/m) 

one sees that 

|zt| < (7m'*(l + |logm|)° < Cm'^, 

which in particular gives us pointwise convergence. As for uniform Holder bounds, take 
s < i so that 

As before, the integral term is bounded by C{t — s)". For the other term, note that 

't-s 



I)Zs 



< C 

< c 



m 



t - s 



m 



A 1 U. 



/3 



< Cit-sf, 
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where we have used the previous point wise estimate on z in the last inequahty. This 
proves that the paths z are uniformly /3-Holder continuous and finishes the proof. □ 

The above proposition shows, for instance, that replacing Brownian motion in our main 
theorem by fractional Brownian motion with Hurst parameter H > 1/2 will not allow for a 
similar statement with non-trivial stochastic area correction. (It is recalled that fractional 
Brownian motion has Holder continuous sample path with exponent arbitrarily close to 
H.) Let us, finally and briefly, discuss a similar statement when Brownian motion in our 
main theorem is replaced by a finite energy path; that is, a path 

7 : [0,r] ^M". 

which may be written as indefinite integral of some function in ([0, T] ,M"), which we 
shall call 7. By Cauchy-Schwarz, such finite energy paths are guaranteed to be 1/2- Holder 
but, in general, one docs not have better Holder regularity. In particular, since the area 
is not continuous in l/2-H61der topology, the above proposition just about fails to cover 
finite energy paths. A direct argument, however, is not difficult. As in the above proof, we 
set 

z := 7o,. - Mxo,. 

and note from the previous argument \zt\ < Cm-*^/^ (1 -|- | log m|)''^^^, uniformly over t G 
[0,r]. We then write 

M 

zt = zt + Jt 

m, 

and take the scalar product in M" with zt, following by integration over [0,T] to see that 

/ \zt\^dt = 7;\zTf+ / {%^)dt 

Jo ^ Jo 

This implies a uniform (in m) L^-bound on i. This of course implies a uniform L^- 
bound on i and thus a uniform 1-variation bound on z. Knowning that z converges to 
zero uniformly on [0,T], it now follows from interpolation that this convergences also takes 
place in p- variation, for any p > 1. Now, the area is a continuous function of the underlying 
paths in p-variation as along as p < 2 and so we can conclude: replacing Brownian motion 
in our main theorem by a finite energy (also known as Cameron-Martin) path will not 
allow for a similar statement with non-trivial stochastic area correction. 

4. Appendix: elements of rough path theory 

A rough path on an interval [0, T] with values in a Banach space V consists of a contin- 
uous function X : [0, T] — )• V, as well as a continuous "second order process" X: [0, T]^ — >■ 
V 'S>V, subject to certain (i) algebraic and (ii) analytic conditions. Towards (i), the be- 
haviour of iterated integrals of smooth paths suggests to impose the algebraic relation 
{"Chen's relation"), 
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assumed to hold for every triple of times {s,u,t). Since Xt^t = 0, it immediately follows 
(take s = u = t) that we also have X^^^ = for every t. One should think of X as postulating 
the value of the quantity 



/ 

J s 



Xs,r dXr — XXg^t ) 



where wc take the right hand side as a definition for the left hand side. We insist that 
knowledge of the path 1 1-^ (Xo,t,Xo,t) already determines the entire second order process 
X. In this sense {X, X) is indeed a path, and not some two-parameter object. 

Note that the algebraic relations are by themselves not sufficient to determine X as a 
function of X. Indeed, for any V <S>V-V8lucd function F, the substitution X^^t i— ^s,t + Ff — 
Fg leaves the left hand side of the above algebraic relation invariant. We will see later on 
how one should interpret such a substitution. The aim of imposing these algebraic relations 
is to ensure that X does indeed have the basic addivity properties of any (reasonable) 
integral when considering it over two adjacent intervals. 

For a G (^, one can denote by P"([0, T],l^) the space of those rough paths (X,X) 
such that 

llx^ll \^s,t\ ^ ii^ii |Xs,t| 

\\X\\a = sup — < OO , ||X||2a = SUp '-pr- < GO . 

s^te[o,T] I* - sr s^telo,T] \t-s\^'^ 

If one ignores the non-linear, algebraic constraint there is a natural way to think of {X, X) 
as element in the Banach space of such maps with (semi-)norm \\X\\a + ||X||2q. However, 
due to the non-linear algebraic relation is not a linear space, although a closed subset 
of the aforementioned Banach space. 

Definition 1. Given rough paths X,Y G X'"([0, T], F), we define the (inhomogenous) 
a-Holder rough path metric 

/~v -<r\ \-^s,t — Ys,t\ , |Xs ( — Ys t\ 

Pa(X,Y):= sup — f— + sup 



s^te[o,T] \t-sr s^te[o,T] \t-s\^" 

Let us note that P"([0, T], V) so becomes a complete, metric space. The perhaps cheap- 
est way to show convergence with respect to this rough path metric is based on interpola- 
tion: in essence, it is enough to establish pointwise convergence in conjunction with uniform 
"rough path" a-H61der bounds. We conclude this part with two important remarks. First, 
we can ask ourselves up to which point the algebraic relations are already sufficient to 
determine X. Assume that we can associate to a given function X two different second 
order processes X and X, and set Gs,t = '^s,t — ^s,t- It then follows immediately that 

Gs,t = Gu,t + Gs,u 

so that in particular Gs.t = Go^t — Gq^s- We conclude that X is in general determined only 
up to the increments of some function F with values in [V V) and Holder continuous 
with exponent 2a. The choice of F does usually matter and there is in general no obvious 
canonical choice. 

The second remark is that this construction can possibly be useful only if a < ^. Indeed, 
if a > ^, then a canonical choice of X is given in terms of the Young integral. Furthermore, 
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it is clear in this case that X must be unique, since any additional increment should be 
2a;-Holder continuous, which is of course only possible if a < ^. This is however not to say 
that X is uniquely determined by X if the latter is smooth, when interpreted as an element 
of P". Indeed, if a < i, F is any 2a-Holder continuous function with values m. V ®V 
and ^s,t = Ft — Fg, then the path (0,X) is a perfectly "legal" element of P", even though 
one cannot get any smoother than the function 0. The impact of perturbing X by some 
F G C^" in the context of differential equations and integration is dramatic: additional 
drift terms in (Lie-bracket) directions can appear; the famous Ito-Stratonovich correction 
is also understood from this picture. The reader may find (much) more in [19l|9] and [8]. 
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